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Motivation 
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Experiments 
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Selection of DDB Model 

• The DDB Model was proposed by Ibrahim, Yang and Przekwas (1993) 

• The DDB is based on the type of droplet deformation observed in 
droplet breakup studies near the leading edge of an airfoil 

- “the liquid droplet is deformed from an initial spherical shape of radius R into an oblate 
spheroid of an ellipsoidal cross section with major semi-axis a and minor semi-axis b" 


• Model governed by a second order ODE 

- Well tried numerical schemes available to do the integration 


• Model prediction of displacement of the center of mass can be 
compared to experimental results 
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DOB Model 



dw 
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Assumptions 

No exchange of heat with surroundings 

Only forces involved: pressure, viscous 
and surface tension 
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DDB Model 
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DDB Model 

Kinetic Energy Term 
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DOB Model 

Potential Energy Term 
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Approximation to the surface area 
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DOB Model 

Potential Energy Term 
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DDB Model 

Work Done by the Pressure 


Approximation to the 
projected area 


Approximation to the 
average pressure 
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DDB Model 


dE _ dW 
dt dt 
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DDB Model 
Pressure Force Term 
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DDB Model 
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viscous force 


surface tension force 


Non-dimensionalization of the equation: 
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Initial Value Problem 
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Input Parameters 
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Parameters for DDB Mode l 
Density of air; viscosity of air 

Density of water; viscosity of water 

Surface tension of water 

Diameter of the droplet 

Slip Velocity 


Parameters for the Numerical Solutio n 
Number of first order ODE 

Initial and final value of the 
independent variable in the interval 
where the solution is evaluated 

Number of steps 

Step size 

Error Tolerance 
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perimental Data 
Rotating Arm 
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Experimental Data 
Data Analysis 

• MATLAB with Digital Imaging Processing tool box was used 
in the data analysis 

• Droplet movement in the horizontal and vertical directions 
tracked frame by frame 

• Program tracks width and height of droplet 

- Ellipse superimposed on the deformed droplet 

- Major and minor semi-axis of the superimposed ellipse 

• Knowing the major semi-axis allows calculation of the center 
of mass vertical displacement for half-droplet 
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Experimental Data 

droplet radius = 516 pm, airfoil velocity = 90 m/sec, airfoil chord = 0,710 m, 
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Experimental Data 

droplet radius = 516 pm, airfoil velocity = 90 m/sec, airfoil chord = 0.710 m, 
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Approach 


• Assume droplet is in quasi-steady equilibrium at each location along the trajectory 

• Solve the model at each location 

- Use experimentally measured slip velocity at each location as model input 

• Compare model prediction to experimental data 
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Results 


Icing Branch 


Center of Mass Oscillation 

location 100, distance from the leading edge of the airfoil = -120.5 millimeters 
droplet radius = 516 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Results 


Center of Mass Oscillation 

location 170, distance from the leading edge of the airfoil = -37.3 millimeters 
droplet radius = 516 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Results 


Icing Branch 


Model Prediction Compared to Experimental Data 
droplet radius = 199 urn, airfoil chord = 0.710 m, airfoil velocity = 50 m/sec. 
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Results 


Model Prediction Compared to Experimental Data 
droplet radius = 287 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Results 


Model Prediction Compared to Experimental Data 
droplet radius = 439 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Results 


Model Prediction Compared to Experimental Data 
droplet radius = 516 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Results 


Model Prediction Compared to Experimental Data 
droplet radius = 685 urn, airfoil chord = 0.710 m, airfoil velocity = 90 m/sec. 
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Conclusions 


• For small and medium size droplets (radius between 200 and 500 pm) the 
model prediction agrees with experimental data. 

• For large droplets (radius larger than 500 pm) the model over-predicts 
displacement of the center of mass by a large margin 

• The increasing deviation between model prediction and experimental data as 
droplet size increases indicates that one or more model assumptions are 
invalid for large droplet sizes 

• The quasi-steady assumption needed to apply the DDB model works well for 
small and medium size droplets 

• The model can be used in the analysis of deformation of small and medium 
droplets from previous experiments 


Icing Branch 


Page 28 


NASA Glenn Research Center 


OF PRESI 


Icing Branch 


Page 29 


NTATION 


NASA Glenn Research Center 


